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Abstract. The explosion of a core collapse supernova drives a powerful shock front 
into the wind from the progenitor star. A layer of shocked circumstellar gas and ejecta 
develops that is subject to hydrodynamic instabilities. The hot gas can be observed 
directly by its X-ray emission, some of which is absorbed and re-radiated at lower 
frequencies by the ejecta and the circumstellar gas. Synchrotron radiation from rel- 
ativistic electrons accelerated at the shock fronts provides information on the mass 
loss density if free-free absorption dominates at early times or the size of the emitting 
region if synchrotron self-absorption dominates. Analysis of the interaction leads to 
information on the density and structure of the ejecta and the circumstellar medium, 
and the abundances in these media. The emphasis here is on the physical processes 
related to the interaction. 

1 Introduction 

The collision of supernova ejecta with dense surrounding gas can generate high 
pressure regions comparable to those found in the emission regions of active 
galactic nuclei. The shock wave interactions are analogous to those occurring 
in 7-ray burst afterglows, although the supernova case involves nonrelativistic 
velocities. Multiwavelength observations, when combined with models for the 
emission, give detailed information on the outer supernova structure and the 
structure of the surrounding mass loss region. The stellar mass loss history 
leading up to the supernova explosion can be deduced from the surrounding 
structure. 

The study of circumstellar interaction around supernovae has benefitted from 
several developments in the observational investigations of supernovae. One is 
the ability to observe supernovae in wavelength regions other than optical. Radio, 
X-ray, infrared, and ultraviolet observations have all been crucial for revealing 
various aspects of circumstellar interaction. Another is the ability to follow the 
optical emission from supernovae to late times with the new generation of large 
telescopes. The emission from circumstellar interaction typically has a lower de- 
cay rate than does the inner emission from gas heated by the initial explosion 
and by radioactivity, so that it eventually comes to dominate the emission. Fi- 
nally, the interest in supernovae generated by SN 1987A and other events has 
increased the rate of supernova discovery, including events which show signs of 
circumstellar interaction from a young age. 
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Many of the observational manifestations of circumstellar interaction are 
treated in other chapters of this volume. The emphasis here is on providing 
the theoretical basis to model the observed phenomena. The structure of the 
freely expanding supernova ejecta and of the surrounding medium provide the 
initial conditions for the interaction and are discussed in the first two sections. 
Section 4 deals with the hydrodynamics of the interaction and section 5 treats 
the emission from hot gas. The effects of supernova radiation on the circum- 
stellar gas are discussed in section 6 and relativistic particles in section 7. The 
discussion and conclusions are in section 8. 

2 Ejecta Structure 

The density structure of a supernova is set up during the first days after the 
explosion. Over this timescale, the pressure forces resulting from the initial ex- 
plosion, and some later power input from radioactivity, become too small to 
change the supernova density structure. The velocity profile tends toward that 
expected for free expansion, v = r/t, and the density of a gas element drops 
as t~ 3 . The pressure drops to a sufficiently low value that it is not a factor in 
considering shock waves propagating in the ejecta as a result of later interaction. 

The density structure of the ejecta is a more complex problem. For core col- 
lapse supernovae, where all the explosion energy is generated at the center of 
the star, the explosion physics lead to an outer, steep power law density profile 
along with an inner region with a relatively flat profile ]l7| , p^ . The outer profile 
is produced by the acceleration of the supernova shock front through the outer 
stellar layers with a rapidly decreasing density. This part of the shock propa- 
gation does not depend on the behavior in the lower layers and the limiting 
structure is described by a self-similar solution. The structure depends on the 
initial structure of the star and thus depends on whether the progenitor star has 
a radiative or a convective envelope In the radiative case, which applies to 
Wolf-Rayet stars and to progenitors like that of SN 1987A, the limiting profile is 
p ^ r -io.2^ j n ^-jje convective case, which applies to red supergiant progenitors, 
the limiting profile is p cx r~ n ' 7 . These are the limiting profiles and the density 
profile over a considerable part of the supernova might be described by a some- 
what flatter profile. Numerical calculations of the explosion of SN 1987A have 
indicated p oc r~( 8 ~ 9 ) in the outer parts of the supernova Q. 

The calculations of the density profile assume an adiabatic flow. As the shock 
wave approaches the surface of the star, radiative diffusion becomes a factor and 
the shock acceleration stops. In the case of SN 1987A, theory as well as observa- 
tions suggest that diffusion is not important until a velocity of ~ 30, 000 km s _1 
is reached Q . A comparable limit applies to the relatively compact progenitors 
of the Wolf-Rayet star explosions Q . In the case of red supergiant progenitors, 
which is the case applicable to most Type II supernovae, a lower velocity limit is 
expected. At the point in the shock evolution where radiative diffusion becomes 
important, the loss of the radiative energy can lead to a large compression of the 
gas, to the point where the thermal pressure becomes important W . There is the 
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possibility of a dense shell at the maximum velocity generated by the supernova. 
The shell may be broken up by Rayleigh- Taylor instabilities [[35| and the final 
outcome of these processes is not known. 

The overall result of these considerations is that the outer part of a core 
collapse supernova can be approximated by a steep power law density profile, or 
/9 e j oc r~ n where n is a constant. After the first few days the outer parts of the 
ejecta expand with constant velocity, V(m) oc r for each mass element, m, so 
that r(m) = V(m)t and p(m) = p (m)(t /t) 3 . Therefore 

Pcj = Po(t/t o y 3 (V t/r) n . (1) 

This expression takes into account the free expansion of the gas. 

Type la supernovae, which are believed to be the thermonuclear explosions 
of white dwarfs, are different because the supernova energy is not all centrally 
generated. A burning front spreads through the star and the energy released 
in this way gives rise to an explosion. The results of a number of explosion 
simulations show that the post-explosion density profile can be approximated 
by an exponential in velocity p7| . 

3 Stellar Mass Loss 

Type II supernovae, with H lines in their spectra, are thought to be the ex- 
plosions of massive stars that have reached the ends of their lives with their 
H envelopes. In most cases, the stars explode as red supergiants, which are 
known to have slow, dense winds. Typical parameters are a mass loss rate of 
M = 10~ 6 — 10~ 4 Mq yr _1 and a wind velocity u w = 5 — 25 km s _1 . If the 
mass loss parameters stay approximately constant leading up to the explosion, 
the circumstellar density is given by 

Pcs = M/(4iru w r 2 ). (2) 

In the late stellar evolutionary phases, the evolution of the stellar core occurs on 
a rapid timescale, but the stellar envelope has a relatively long dynamical time, 
which can stabilize the mass loss properties. 

The mechanisms by which mass is lost in the red supergiant phase are poorly 
understood, but some insight into the wind properties can be gained by consid- 
ering observations of the winds. VY CMa, thought to have a zero-age main 
sequence mass of 30 — 40 Mq ^J, a mass loss rate ~ 3 x 10~ 4 Mq yr -1 , 
and u w « 39 km s _1 p5| , |69| , is an especially well-observed case. HST imag- 
ing shows that the density profile is approximately r~ 2 over the radius range 
3 x 10 16 — 1.4 x 10 17 cm but that there is considerable structure superposed on 
this profile, including knots and filaments p9f . Observations of masers imply the 
presence of clumps with densities up to ~ 5 x 10 9 cm -3 and suggest the pres- 
ence of an expanding disk seen at an oblique angle f37f . VY CMa is an extreme 
mass loss object, but there is evidence for irregular mass loss in a Orionis and 
other red supergiants |75|,f76|. Among lower mass stars, the bipolar structure of 
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planetary nebulae has been widely attributed to a high equatorial density in the 
red giant wind. In one case where we can observe the circumstellar surroundings 
of a supernova, SN 1987A, an axisymmetric structure similar to that found in 
planetary nebulae has been observed [||. 

Some massive stars, either through individual mass loss or through binary 
interaction, lose their H envelopes entirely and become Wolf-Rayet stars. Their 
ultimate explosions are thought to be observed as Type lb and Ic supernovae. 
The progenitor stars are more compact in this case and have faster winds, with 
u w = 1, 000 - 2, 500 km s" 1 and M = 10~ 6 - 10~ 4 M yr^ 1 fcj. The fast wind 
can create a bubble in the surrounding medium, which is typically the slow wind 
from a previous evolutionary phase, and the resulting shells have been observed 
around a number of Wolf-Rayet stars |45| . Their typical radii are a few pc. 

The progenitors of Type la supernovae are not known, so that observations 
of circumstellar interaction as well as other signatures of the companion star po- 
tentially contain important clues on the progenitor question ||,U,^9| . In models 
where a white dwarf accretes mass from a companion star, the wind from the 
companion can provide a relatively dense circumstellar medium. In the case of a 
double degenerate progenitor, there may be a disk around the coalesced object, 
but the interaction is with the interstellar medium on larger scales. There have 
been no detections of circumstellar interaction around Type la supernovae, so 
we do not discuss them further in this review, although the physical processes 
described here should apply if any surrounding wind is present. 

4 Hydrodynamics 

When the radiation dominated shock front in a supernova nears the stellar sur- 
face, a radiative precursor to the shock forms when the radiative diffusion time 
is comparable to the propagation time. There is radiative acceleration of the gas 
and the shock disappears when optical depth ~ unity is reached Q. The fact 
that the velocity decreases with radius implies that the shock will re-form as 
a viscous shock in the circumstellar wind. This occurs when the supernova has 
approximately doubled in radius. 

The interaction of the ejecta, expanding with velocity > 10 4 km s _1 , and 
the nearly stationary circumstellar medium results in a reverse shock wave prop- 
agating inwards (in mass), and an outgoing circumstellar shock. The density in 
the circumstellar gas is given by equation (Q) . As discussed above, hydrodynam- 
ical calculations show that to a good approximation the ejecta density can be 
described by equation ([j]). A useful similarity solution for the interaction can 
then be found flUQ. Here we sketch a simple derivation. More details can be 
found in these papers, as well as in the review 

Assume that the shocked gas can be treated as a thin shell with mass M s , 
velocity V s , and radius R s . Balancing the ram pressure from the circumstellar 
gas and the impacting ejecta, the momentum equation for the shocked shell of 
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circumstellar gas and ejecta is 



M s ^ = a-*r 2 Mv - V a f - Pcs V s 2 }. 



(3) 



Here, M s is the sum of the mass of the shocked ejecta and circumstellar gas. 
The swept up mass behind the circumstellar shock is M cs — AIR s /u w , and 
that behind the reverse shock M rev = Air t^V™ (t / i? s )™~ 3 / (n — 3), assuming that 
Rs >> R p , the radius of the progenitor. With V — R s /t we obtain 



M , 4tt Po tl V a n t"- 3 
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This equation has the power law solution 
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(5) 



The form of this similarity solution can be written down directly by dimensional 
analysis from the only two independent quantities available, p t 3 , V" and M /u w . 
The solution applies after a few expansion times, when the initial radius has been 
'forgotten.' The requirement of a finite energy in the flow implies n > 5. More 
accurate similarity solutions, taking the structure within the shell into account, 
are given ha |6 64 . In general, these solutions differ by less than ~ 30% from the 
thin shell approximation. 

The maximum ejecta velocity close to the reverse shock depends on time as 
V = R s /t oc £ -1 /( Tl-2 ). The velocity of the circumstellar shock, dR s /dt, in terms 
of V is V s = V(n — 3)/(n — 2) and the reverse shock velocity, V ICV = V — V s — 
V/(n — 2). Assuming cosmic abundances and equipartition between ions and 
electrons, the temperature of the shocked circumstellar gas is 



1.36 x 10 a 



n — 3 
n-2 



V 



10 4 km s" 



K 



and at the reverse shock 



(n-3) 2 ' 

The time scale for equipartition between electrons and ions is 

1.5 



10 7 cm- 



(6) 



(7) 



(8) 



One finds that the reverse shock is marginally in equipartition, unless the tem- 
perature is > 5 x 10 s K. The ion temperature behind the circumstellar shock is 




Fig. 1. Density and temperature structure of the reverse and circumstellar shocks for 
n — 7 and a velocity of 2.5 x 10 4 km s -1 at 10 days. Both shocks are assumed to be 
adiabatic. Because of the slow Coulomb equipartition the electron temperature (dotted 
line) is much lower than the ion temperature (solid line) behind the circumstellar shock. 



> 6 x 10 9 K for Vi > 1.5, and the density a factor > 4 lower than behind the 
reverse shock. Ion-electron collisions are therefore ineffective, and T c << Ti on , 
unless efficient plasma instabilities heat the electrons collisionlessly (fig. [j]). 

For typical parameters, the electron temperatures of the two shocks are very 
different, ~ (1 — 3) x 10 9 K for the circumstellar shock and 10 7 — 5 x 10 s K for the 
reverse shock, depending on n. The radiation from the reverse shock is mainly 
below ~ 20 keV, while that from the circumstellar shock is above ~ 50 keV. 

The density behind the reverse shock is 

_ (n-4)(n-3) 

Prcv — ^ Pes 

and is much higher than behind the circumstellar shock for n > 7. There is 
a drop in density across the contact discontinuity, moving from the shocked 
ejecta to the circumstellar medium (see fig. ||). The fact that low density gas is 
decelerating higher density gas leads to a Rayleigh- Taylor instability. Chevalier, 
Blondin & Emmering [ |l8| have calculated the structure using a two-dimensional 
PPM (piecewise parabolic method) hydrodynamic code. They indeed find that 
instabilities develop, with dense, shocked ejecta gas penetrating into the hotter, 
low density shocked circumstellar gas (fig. ^). The instability mainly distorts 
the contact surface, and does not seriously affect the general dynamics. The 
calculation assumed that cooling is not important. If the gas at the reverse shock 
cools efficiently, the extent of the instability is similar, although the Raylcigh- 
Taylor fingers are narrower Jl3| . 

In view of the evidence for dense equatorial winds from red supergiant stars, 
Blondin, Lundqvist, & Chevalier J| simulated the interaction of a supernova with 
such a wind. They found that for relatively small values of the angular density 



(9) 
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Fig. 2. Two-dimensional calculation of the shock structure for a supernova with n = 6 
in a stellar wind (courtesy John Blondin). 



gradient, the asymmetry in the interaction shell is greater than, but close to, 
that expected from purely radial motion. If there is an especially low density 
close to the pole, the flow qualitatively changes and a protrusion emerges along 
the axis, extending to 2 — 4 times the radius of the main shell. Protrusions have 
been observed in the probable supernova remnant 41.9+58 in M82, although the 
nature of the explosion is not clear in this case Q . 

In addition to asymmetric winds, there is evidence for supernova shock waves 
interacting with clumps of gas in the wind, as have been observed in some red 
supergiant winds (section ||) . In some cases the clumps can be observed by their 
very narrow lines in supernova spectra, as in Type Iln supernovae. The velocity of 
a shock wave driven into a clump, v c , can be estimated by approximate pressure 
balance v c ~ v s (p s / Pc) 1 ^ 2 , where v s is the shock velocity in the smooth wind 
with density p s and p c is the clump density. The lower shock velocity and higher 
density can lead to radiative cooling of the clump shock although the main shock 
is non-radiative. Optical line emission of intermediate velocity observed in Type 
Iln (narrow line) supernovae like SN 1978K, SN 1988Z, and SN 1995N can be 
explained in this way |2l] , |22] , ^4]| . 

The presence of many clumps can affect the hydrodynamics of the interaction. 
Jun, Jones, & Norman |49|j found that propagation in a region with clumps gives 
rise to widespread turbulence in the shocked region between the forward shock 
and the reverse shock, whereas the turbulence is confined to a region near the 
reverse shock for the non-clump case (fig. |^). Their simulations are for interaction 
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with a constant density medium, but the same probably holds true for interaction 
with a circumstcllar wind. 



5 Emission from the hot gas 

During the first month, radiation from the supernova photosphere is strong 
enough for Compton scattering to be the main cooling process for the circum- 
stellar shock. The photosphcric photons have energies ~ 3 kT e ff 1 — 10 eV. 
The optical depth to electron scattering behind the circumstellar shock is 

r c = 0.18 M_ s u-{ Vf 1 t~l ys . (10) 

A fraction of the photospheric photons will scatter N times in the hot 
gas. In each scattering the photon increases its energy by a factor Avjv s» 
4 kT c /m c c 2 > 1- The multiple scattering creates a power law continuum that 
may reach as far up in energy as the X-ray regime. If relativistic effects can be 
ignored (T c < 10 9 K), the spectral index is |37) 

° "^' 2 ln[5(0.9228 -lnr o )]l 1/2 -|. (11) 



kT c 2 wjj 2 

Typically, 1 < a < 3. This type of emission may have been observed in the ultra- 



violet emission from SN 1979C 65 37|. For T c > 10 9 K relativistic effects become 
important and considerably increase the cooling |5^] . 

One can estimate the free-free luminosity from the circumstellar and reverse 
shocks from 

U = 4tt / A s (T c )n 2 c r 2 dr « A B (TJ . MiPi .. . (12) 

where the index i refers to quantities connected either with the reverse shock 
or circumstellar shock. The density behind the circumstellar shock is p cs = 
4 po = M/(iru w R 2 ). The swept up mass behind the circumstellar shock is 
M cs = MR s /u w and that behind the reverse shock M rcv — (n — 4)M cs /2. With 
A s = 2.4 x W- 27 g s T c 5 , we get 

fc-MxlB-fcG.^)'^)" 1 = g8 -, (13) 

where gg is the free-free Gaunt factor, including relativistic effects. For the re- 
verse shock C n = (n — 3)(n — 4) 2 /4(n — 2), and for the circumstellar shock 
C n = 1. This assumes electron-ion equipartition, which is highly questionable 
for the circumstellar shock (see fig. Q). Because of occultation by the ejecta only 
half of the above luminosity escapes outward. 

At T c < 2 x 10 7 K, line emission increases the cooling rate and A w 3.4 x 
10 -23 T e ^°' 67 erg s _1 cm 3 . If the temperature of the reverse shock falls below 
~ 2 x 10 7 K, a thermal instability may occur and the gas cools to < 10 4 K, where 
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photoelectric heating from the shocks balances the cooling. Using t coo \ — 3kT e /A, 
one obtains for the cooling time 

605 / v ci \ 5 34 (m-^Y 1 ( t V 

CODl ~ (n - 3)(n - 4)(n - 2) 3 - 34 ^ 10 4 km s" 1 J ^ Uu)1 J ^days J yS ' 

(14) 

assuming solar abundances [[id] . From this expression it is clear that the cooling 
time is very sensitive to the density gradient, as well as the shock velocity and 
mass loss rate. SNe with high mass loss rates, like SN 1993J, generally have 
radiative reverse shocks for > 100 days, while SNe with lower mass loss rates, 
like the Type IIP SN 1999em, have adiabatic shocks from early times. 

The most important effect of the cooling is that the cool gas may absorb 
most of the emission from the reverse shock. Therefore, in spite of the higher 
intrinsic luminosity of the reverse shock, little of this will be directly observable. 
The column density of the cool gas is given by iV coo i = Af rcv /(47r i? 2 m p ), or 

(£) (tom^t)" (look)" ™-' 2 ' 

(15) 

Because the threshold energy due to photoelectric absorption is related to iV C ooi 
by E(t = 1) = 1.2(iV coo i/10 22 cm -2 ) 3 / 8 keV, it is clear that the emission from 
the reverse shock is strongly affected by the cool shell, and a transition from 
optically thick to optically thin is expected during the first months, or year. As 
an illustration, we show in figure || the calculated X-ray spectrum at 10 days and 



at 200 days for SN 1993 J 41 .At early epochs the spectrum is dominated by the 
very hard spectrum from the circumstellar shock, which reaches out to > 100 
keV. At later epochs the soft spectrum from the reverse shock penetrates the 
cool shell, and the line dominated emission from the cooling gas dominates. 
If cooling, the total energy emitted from the reverse shock is 

T ^tf l nV* (n-3)(n-i) MV* 
L ICV - 4irK s - p cj V rcv = 



2 I C3 rov 4(n-2) 3 



= 1.6 x 10 41 ^(^s - M- 5 u-±V? erg b" 1 . (16) 

For high M/u w the luminosity from the reverse shock may contribute apprecia- 
bly, or even dominate, the bolometric luminosity. 

Because V oc t" 1 ^ n ^ 2 \ L rcv oc £~ 3 /("- 2 ) in the cooling case. Although the 
total luminosity is likely to decrease in the cooling case, the increasing trans- 
parency of the cool shell, r coo i oc can cause the observed flux in energy 
bands close to the low energy cutoff, E(t = 1), to increase with time, as was 
seen, e.g., in SN 1993J (figure |). 

Because of the low temperature the spectrum of the reverse shock is dom- 
inated by line emission from metals (fig. ||). An important point is that the 
observed spectrum is formed in gas with widely different temperatures, varying 
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Fig. 3. X-ray spectrum of SN 1993J at 10 days and at 200 days. At 10 days the free- free 
emission from the outer shock dominates, while at 200 days the cool shell is transparent 
enough for the line dominated spectrum from the reverse shock to dominate instead. 
Right: The solid lines give the luminosity in the 0.1-2.4, 1-10, and 50-100 keV bands, 
corrected for absorption, as a function of time, while the dotted lines give the total 
emitted luminosity from the reverse and circumstellar shocks jll| 

from the reverse shock temperature to ~ 10 4 K. A spectral analysis based on a 
one temperature model can be misleading. 

Chugai |l9| has proposed that the X-ray emission from the Type Iln SN 
1986J is the result of the forward shock front moving into clumps, as opposed to 
the reverse shock emission. One way to distinguish these cases is by the width 
of line emission; emission from the reverse shock wave is expected to be broad. 
It has not yet been possible to carry through with this test |l6| . 

Another way of observing the hot gas is through the emission from collision- 
ally heated dust grains in the gas p0| ] . Dust formation in the rapidly expanding 
ejecta is unlikely, so the forward shock front must be considered. Evaporation 
by the supernova radiation creates a dust-free zone around the supernova (see 
section 6.3). The time for the supernova shock wave to reach the dust depends on 
the supernova luminosity and the shock velocity; it is probably at least several 
years. The infrared luminosity from dust can be up to ~ 100 times the X-ray 
luminosity of the hot gas for typical parameters, and the dust temperature is a 
measure of the density of the gas . If the X-ray emission from a supernova 
like SN 1986J is from circumstellar clumps that are out in the region where dust 
is present, there is the possibility of a large infrared luminosity. 

6 Radiative heating and re-emission 
6.1 Soft X-ray burst and circumstellar gas 

The earliest form of circumstellar interaction occurs at shock break-out. As the 
shock approaches the surface, radiation leaks out on a time scale of less than an 
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hour. The color temperature of the radiation is ~ (1 — 5) x 10 5 K and the energy 
(1 - 10) x 10 48 ergs §i]@||,|§. This burSt ° f EUV ( extreme ultraviolet) 



and soft X-rays ionizes and heats the circumstellar medium on a time scale of 
a few hours. In addition, the momentum of the radiation may accelerate the 
circumstellar gas to a high velocity. Most of the emission at energies > 100 eV 
is emitted during the first few hours, and after 24 hours little ionizing energy 
remains. 

The radiative effects of the soft X-ray burst were most clearly seen from 
the ring of SN 1987A, where a number of narrow emission lines from highly 
ionized species, like N III-N V, were first seen in the U V |4^ , [72| . Later, a forest 
of lines came to dominate also the optical spectrum ]78|. Imaging with HST 
(e.g., f^jH) showed that the lines originated in the now famous circumstellar 
ring of SN 1987A at a distance of ~ 200 light days from the SN. The presence of 
highly ionized gas implied that the gas must have been ionized and heated by the 
radiation at shock break-out. Because of the finite light travel time across the 
ring, the observed total emission from the ring is a convolution of the emission 
at different epochs from the various part of the ring. Detailed modeling ]55| . [5G[ 
shows that while the ionization of the ring occurs on the time scale of the soft 
X-ray burst, the gas recombines and cools on a time scale of years, explaining 
the persistence of the emission decades after the explosion. The observed line 
emission provides sensitive diagnostics of both the properties of the soft X-ray 
burst, and the density, temperature and abundances of the gas in the ring. 
In particular, the radiation temperature must have reached ~ 10 6 K, in good 
agreement with the most detailed recent modeling of the shock break-out ||. 
Narrow emission lines are not unique to SN 1987A, but have also been observed 
for several other SNe, in particular several Type Iln SNe, such as SN 1995N Ej 
and SN 1998S §§. 

The soft X-ray burst may also pre-accelerate the gas in front of the shock. In 
the conservative case that Thompson scattering dominates, the gas immediately 
in front of the shock will be accelerated to 

T/=1.4xl0 3 (—£ ^ ( ^ T ) 2 (—) 2 kms- 1 , (17) 

V. 10 48 ergs J \1 x 10 4 km s" 1 J V d ays / V ' 

where E is the total radiative energy in the burst. Line absorption may fur- 
ther boost this p9[| . If the gas is pre-accelerated, the line widths are expected 
to decrease with time. After about one expansion time (~ R p /V) the reverse 
and circumstellar shocks are fully developed, and the radiation from these will 
dominate the properties of the circumstellar gas. The fraction of this emission 
going inward is absorbed by the ejecta and there re-emitted as optical and UV 
radiation ]37| , |38| |. 

The X-ray emission from the shocks ionizes and heats both the circumstel- 
lar medium and the SN ejecta. Observationally, these components are distin- 
guished easily by the different velocities. The circumstellar component is ex- 
pected to have velocities typical of the progenitor winds, i.e., < 1000 km s _1 , 
while the ejecta have considerably higher velocities. The density is likely to 
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be of the order of the wind density 10 5 — 10 cm" 3 , or higher if clumping is 
important. The ionizing X-ray flux depends strongly on how much of the flux 
from the reverse shock can penetrate the cool shell. The state of ionization in 
the circumstellar gas is characterized by the value of the ionization parameter, 
C = L cs /(r 2 n) = 10 2 (L CS /10 40 erg s^fr/lO 16 cm)~ 2 (n/10 6 cm" 3 )- 1 @. The 
comparatively high value of £ « 10 — 10 3 explains the presence of narrow coronal 



lines of [Fe V-XI] seen in objects like SN 1995N |44 



6.2 SN ejecta 

The ingoing X-ray flux from the reverse shock ionizes the outer parts of the 
ejecta. The state of highest ionization therefore is close to the shock, with a 
gradually lower degree of ionization inwards. Unless clumping in the ejecta is 
important, the ejecta density is ~ 10 6 — 10 8 cm -3 . In the left panel of figure 
|] we show temperature and ionization structure of the ejecta, as well as the 
emissivity of the most important lines. The temperature close to the shock is 
~ 3 x 10 4 K. Calculations show that most of the emission here is emitted as U V 
lines of highly ionized ions, like Lya, C IIITV, N III-V, and O III- VI. Inside the 
ionized shell there is an extended partially ionized zone, similar to that present 
in the broad emission line regions of AGN's. Most of the emission here comes 
from Balmer lines. 

As we have already discussed, the outgoing flux from the reverse shock is to 
a large extent absorbed by the cool shell between reverse shock and the contact 
discontinuity if radiative cooling has been important. The whole region behind 
the reverse shock is in approximate pressure balance, and the density of this 
gas is therefore be a factor ~ 4T rov /T coo i w 10 3 — 10 4 higher than that of the 
ejecta. Because of the high density, the gas is only be partially ionized and the 
temperature only (5 — 8) x 10 3 K. Most of the emission comes out as Balmer 
lines, Mg II and Fe II lines (figure^, right panel). The thickness of the emitting 
region is also very small, ~3x 10 12 cm. In one dimensional models, the velocity 
is marginally smaller than the highest ejecta velocities. Instabilities in the shock 
are, however, likely to erase this difference. 

An important diagnostic of the emission from the cool shell and the ejecta 
is the Ha line. This line arises as a result of recombination and collisional 
excitation. In [^6), it is shown that ~ 1% of the reverse shock luminosity is 
emitted as Ha, fairly independent of density and other parameters. Observations 
of this line permit us to follow the total luminosity from the reverse shock, 
complementary to the X-ray observations. In SN 1993J, the Ha line had the box- 
like shape that is expected for shocked, cooled ejecta |6^^l[. The top of the line 
showed structure that varied with time; this could be related to hydrodynamic 
instabilities of the reverse shocked gas. 



6.3 Interaction with dust 

The winds from red supergiant stars are known to contain dust, so that infrared 
emission from radiatively heated dust and dust scattering of supernova light 
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Fig. 4. Structure of the ejecta and cool shell ionized by the reverse shock at 500 days 
for parameters appropriate to SN 1993J (M = 5 x lO~ 5 Af yr _1 for u = 10 km s _1 ). 
Upper panels show the temperature and ionization of the ejecta (left panels) and the 
cool shell (right panels) , while the lower panels show the corresponding luminosities per 
unit distance. Note the different length scales in the two panels. The ejecta region has 
low density, high temperature and ionization, while the cool shell has a high density, 
is extremely thin, has a low temperature, and is only partially ionized. 



might be expected. The dust temperature is determined by a balance between 
radiative heating and emission from the grain surface. If the supernova radiates 
like a blackbody of temperature T sn and if the dust absorption efficiency varies 
as A -9 , the grain temperature is T„ = T sn W /1 ^ 4+9 \ where W is the dilution 
factor for the supernova radiation p8fl . For 5 = 1, we have 

Tg = 280 (dWi<) (lO«ergss-0 (l0^) K ' (18) 

where L sn is the supernova luminosity. 

Because typical evaporation temperatures for grain materials are in the range 
1000 — 1500 K, the dust near the supernova is evaporated. The radius out to 
which dust is evaporated, r v , is probably determined by the luminosity at the 
time of shock breakout, when L sn can be > 10 44 ergs s _1 . SN 1979C and SN 
1980K were detected as infrared sources, and Dwek |28| estimates r„«3x 10 17 
cm in these cases. The infrared light curve can be calculated from T g , taking into 
account light travel time effects. If the characteristic time that the supernova is 
bright is short compared to 2r„/e, there is a plateau phase until a time 2r v /c. The 
infrared flux then drops as t~ 2 for a highly extended wind and drops more rapidly 
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if the cutoff in the wind, r w , is close to r v . For SN 1979C and SN 1980K, Dwek 
[ p8| found r w < 10 18 cm. The ratio of the infrared emitted energy to the total 
emitted energy gives an estimate of the optical depth through the dusty wind. 
The optical depths for SN 1979C and SN 1980K are ~ 0.3 and 0.03, respectively, 
leading to minimum shell masses of ~ 1 — 5 Mq and ~ 0.1 — 0.4 M for the 
two supernovae Q . The corresponding mass loss rates are consistent with those 
derived from radio observations [^4| . Infrared emission that can be attributed to 
radiatively heated dust has recently been observed from SN 1998S |Q. Infrared 
dust echoes have the potential to give information on dust composition and 
distribution ||^,^|, but the observations do not yet exist to address these issues. 

The dust grains that give an infrared echo can also give a scattered light echo 
[ [To| . The ratio between the total scattered light and the infrared light depends 
on the albedo of the dust. The light curve decreases with time before 2r v /c in 
this case because of the strong forward scattering of typical grains. Chevalier 
[ [To| failed to find good evidence for scattered light echoes from SN 1979C and 
SN 1980K, although such echoes might have been observable. The implication 
may be that the dust grains have small albedos. The problem with detecting 
scattered light echoes is that any echo light may be dominated by light coming 
directly from the circumstellar interaction. Roscherr & Schaefer |}8| examined 
the late emission from two Type Iln supernovae and determined that it was due 
to shock interaction. Because of the ability to spatially resolve the emission in SN 
1987A, it has been possible to observe scattered light from the wind interaction 
nebula p3[ and possibly a shell at the outer extent of the red supergiant wind 




7 Relativistic Particles 



Unambiguous evidence for the presence of relativistic electrons comes from radio 
observations of SNe. A characteristic is the wavelength-dependent turn-on of the 
radio emission ((73j; chapter by Sramek), first seen at short wavelengths, and 
later at longer wavelengths. This behavior is interpreted as a result of decreasing 
absorption due to the expanding emitting region |5). 

Depending on the magnetic field and the density of the circumstellar medium, 
the absorption may be produced either by free-free absorption in the surround- 
ing thermal gas, or by synchrotron self-absorption by the same electrons that 
are responsible for the emission. The relativistic electrons are believed to be 
produced close to the interaction region, which provides an ideal environment 
for the acceleration of relativistic particles. The details of the acceleration and 
injection efficiency are still not well understood (see, e.g., [|6j and references 
therein). Here we just parameterize the injection spectrum with the power law 
index p\ and an efficiency, ij, in terms of the postshock energy density. Without 
radiation or collisional losses the spectral index of the synchrotron emission will 
then be a = (p — l)/2, where flux oc v~ a . Diffusive acceleration predicts that 
Pi = 2 in the test particle limit. If the particle acceleration is very efficient and 
nonlinear effects are important, the electron spectrum can be steeper (e.g., pH). 
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For free-free absorption, the optical depth rg — J_ Ksn e nidr from the radio 
emitting region close to the shock through the circumstellar medium decreases as 
the shock wave expands, explaining the radio turn-on. Assuming a fully ionized 
wind with constant mass loss rate and velocity, so that equation (^) applies, the 
free-free optical depth at wavelength A is 

r s (X) « 7.1 x 10 2 A 2 T-^V 4 -\l (19) 

where M_5 is the mass loss rate in units of 10 -5 Mq yr _1 , u w \ the wind velocity 
in units of 10 km s _1 , and T5 the temperature of the circumstellar gas in 10 5 K. 
From the radio light curve, or spectrum, the epoch of ts = 1 can be estimated 
for a given wavelength, and from the line widths in the optical spectrum the 
maximum expansion velocity, V, can be obtained. Because the effects of the 
radiation from the supernova have to be estimated from models of the circum- 
stellar medium, the temperature in the gas is uncertain. Calculations show that 
initially the radiation heats the gas to T c « 10 5 K |54j| . T c then decreases with 
time, and after a year T « (1.5 — 3) x 10 4 K. In addition, the medium may 
recombine, which further decreases the free-free absorption. From t[r(X)s = 1] 
the ratio M/u w can be calculated. Because M/u w oc Tc^x" 1 , errors in T e and 
x c may lead to large errors in M. If the medium is clumpy, equation ( |l9| ) may 
lead to an overestimate of M /u w . 

Under special circumstances (see below), synchrotron self-absorption (SSA) 
by the same relativistic electrons emitting the synchrotron radiation may be 
important [ [70|Jl^ , ^0] | . The emissivity of the synchrotron plasma is given by j(X) = 
const. X a B~ a N TC \ Q, while the optical depth to self- absorption is given by 

t s = const. A (5/2)+Q B (3/2)+Q ^ rc i . (20) 

Here A^ re ] is the column density of relativistic electrons and B the magnetic 
field. The flux from a disk with radius R s of relativistic electrons is, including 
the effect of SSA, given by 

F u (X)ocR 2 s S(X)[l-e-^% (21) 

where 5(A) = j(X)/n(X) = const. X~ 5 ^ 2 B~ 1 / 2 is the source function. In the opti- 
cally thick limit we therefore have ^(A) w const, i? 2 A -5 / 2 i? -1 / 2 , independent 
of iV re i. A fit of this part of the spectrum therefore gives the quantity i? 2 B^ 1 / 2 . 
The break of the spectrum determines the wavelength of optical depth unity, 
A(r s = 1). Equation ( ^o| ) therefore gives a second condition on iJ 3 / 2 + Q N Te \. If 
R s (t) is known in some independent way, one can determine both the magnetic 
field and the column density of relativistic electrons, independent of assump- 
tions about equipartition, etc. In some cases, most notably for SN 1993J, the 
shock radius, R s , can be determined directly from VLBI observations. If this is 
not possible, an alternative is from observations of the maximum ejecta veloc- 
ity seen in, e.g., the Ha line, which should reflect the velocity of the gas close 
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Fig. 5. Magnetic field (left) and density of relativistic electrons (right) as a function 
of the shock radius for SN 1993 J. The dashed lines show the expected evolution if the 
magnetic energy density and relativistic particle density scale with the thermal energy 
density, B 2 /8tt oc pV 2 oc n TC \ oc t~ 2 , while the dotted lines show the case when B oc r _1 
and n rc i oc r~ 2 Q. 



to the shock. Because the SN expands homologously, R s = V max t. A fit of the 
spectrum at a given epoch can therefore yield both B and N re \ independently. 
From observations at several epochs the evolution of these quantities can then 
be determined. 

Although the injected electron spectrum from the shock is likely to be a 
power law with pi sa 2 (cv rs 0.5), the integrated electron spectrum is affected 
by various loss processes. Most important, the synchrotron cooling time scale of 
an electron with Lorentz factor 7 is t syn w9 x 10 3 7 _1 -B -2 days. This especially 
affects the high energy electrons, steepening the index of the electron spectrum 
by one unit, p = Pi + 1 ~ 3 (a « 1). Inverse Compton losses have a similar effect. 
At low energy, Coulomb losses may be important, causing the electron spectrum 
to flatten. 

The best radio observations of any SN were obtained for SN 1993J. This SN 
was observed from the very beginning until late epochs with the VLA at wave- 
lengths between 1.3-90 cm JT7| , producing a set of beautiful light curves (see 
also Chapter by Sramek). In addition, the SN was observed with VLBI fl57|j5"8| ,|2]] , 
resulting in an impressive sequence of images in which the radio emitting plasma 
could be directly observed. These images showed a remarkable degree of sym- 
metry and clearly resolved the shell of emitting electrons. The evolution of the 
radius of the radio emitting shell could be well fitted by i? s cx i 86 , implying a 
deceleration of the shock front. 

From a fit of the observed spectra for the different epochs the magnetic 
field and number of relativistic electrons could be determined for each epoch, as 
described above p(| . In figure [| we show the evolution of B and particle density 
n Te i, plotted as a function of the shock radius. The most remarkable thing is 
the smooth evolution of these quantities, showing that B w 6.4(i? s /10 16 cm) -1 
G, and n re i oc pV 2 oc t~ 2 , the thermal energy density behind the shock. The 
magnetic field is close to equipartition, B 2 /8ir w 0.14 pV 2 , much higher than 
expected if the circumstellar magnetic field, of the order of a few mG, was just 




Fig. 6. Observed @ and model radio light curves of SN 1993 J @. 

compressed, and strongly argues for field amplification, similar to what has been 
seen in simulations . Contrary to earlier, simplified models for SN 1993J based 
on free-free absorption only jll],(77) , the circumstellar density was consistent with 
p oc r -2 . 

In figure |], we show the excellent fit of the resulting light curves. The high 
values of B implied that synchrotron cooling was important throughout most 
of the evolution for the electrons responsible for the cm emission, and also for 
the 21 cm emission before ~ 100 days. At early epochs, Coulomb losses were 
important for the low energy electrons. The injected electron spectrum was best 
fitted with pi = 2.1. 

The form of the light curves can be understood if, for simplicity, we assume 
equipartition, so that B 2 /8it = rjpV 2 . With p oc (M/u) R~ 2 and V s oc R s /t, we 
find that B oc (M/tt) 1 / 2 ^ 1 . The optically thick part is therefore given by 

F V {X) oc R 2 S \- 5 / 2 B-^ 2 oc (Af/ur^A^/V 5 "- 14 ^"- 2 ), (22) 

since R s oc t("~ 3 )/("~ 2 ). For large n, we get F„(X) oc t 5 ^ 2 . An additional curvature 
of the spectrum is produced by free-free absorption in the wind, although this 
only affects the spectrum at early epochs. 

In the optically thin limit, F V (X) oc R 2 j(X) oc R 2 X a B 1+a N lcl . If losses are 
unimportant, iV re i tot = 47r_R 2 A r rc i, the total number of relativistic electron, may 
either be assumed to be proportional to the total mass, if a fixed fraction of the 
shocked electrons are accelerated, or be proportional to the swept up thermal 
energy. In the first case, A^^ot oc MR s /u w , while in the second iV re i,tot °c 
MR S V 2 /u w , so that in general A^i.tot MR s V 2e /u w , where e = or 1 in 
these two cases. Therefore, F V (X) oc M/u R s V 2e X a B 1+a . If the B-field is in 
equipartition, as above, and using V = (n — 3)/(n — 2)R s /t oc t _1 /(™~ 2 ) we find 

F„(A) OC (M/u w Y 3+a ^ 2 A af-a-(l+2e)/(n-2)_ (23) 
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Fig. 7. Peak luminosity and corresponding epoch for the well-observed radio SNe. The 
dashed lines give curves of constant expansion velocity, assuming SSA (see fH|). 



If synchrotron cooling is important, a similar type of expression can be derived 
[ [Io| . The main thing to note is, however, that the optically thin emission is ex- 
pected to be proportional to the mass loss rate, M/u w , and that the decline 
rate depends on whether the number of relativistic particles scale with the num- 
ber density or the thermal energy of the shocked gas, as well as spectral index. 
Observations of the decline rate can therefore test these possibilities. 

Although a self-consistent model can be developed for SN 1993J and other 
radio supernovae, the modeling of SN 1986J and related objects has been unclear. 
Weiler, Panagia, & Sramek Q proposed a model for SN 1986J in which thermal 
absorbing material is mixed in with the nonthermal emission; one possibility for 
this is a very irregular shocked emission region. The absorption is described 
in a parameterized way and is not related to a quantitative supernova model. 
Chugai & Belous [^0| propose a model in which the absorption is by clumps. The 
narrow line optical emission implies the presence of clumps, but they are different 
from those required for the radio absorption. The possible presence of clumps 
and irregularities introduces uncertainties into models for the radio emission, 
although rough estimates of the circumstellar density can still be obtained. 

The relative importance of SSA and free- free absorption depends on a number 
of parameters and we refer to |||o| for a more detailed discussion. The most 
important of these are the mass loss rate, M/u w , the shock velocity, V s , and 
the circumstellar temperature, T c . In general, a high shock velocity and a high 
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circumstellar temperature favor SSA, while a high mass loss rate favors free-free 
absorption. 

Assuming SSA, an interesting expression for the velocity of the shock can be 
derived which can be tested against the observations. If we assume equipartition, 
iV rc i oc pV 2 R s and B 2 /8ir oc pV 2 we have AT rc i oc B 2 R S . Using this in equation 
( p0| ) we get t s = const \ 5 / 2 + a B 7 / 2 + a j{ s . The peak in the light curve is given 
by t s rs 1. If we approximate the flux at this point by the optically thick expres- 
sion equation (p2) and solve for B we get B oc F v {\)~ 2 i? 4 A -5 . Inserting this 
expression in the condition r s » 1, we find Rl 5+4a Fu(\)~~ 7 ~ 2a \- 15 - ia ~ const. 
With V s = (n - 3)/(n - 2)R s /t we finally have 

V s « const ^(A)( 7 + 2 «)/( 15 + 4 «) A i" 1 (24) 

where all parameters refer to their values at the peak of the light curve. Using 
this expression, we can plot lines of constant shock velocity into a diagram with 
peak radio luminosity versus time of peak flux, assuming that SSA dominates 
(figure ^). The positions of the lines depend only weakly on the equipartition 
assumption. Each SN can now be placed in this diagram to give a predicted 
shock velocity. If this is lower than the observed value (as measured by VLBI 
or from line profiles) SSA gives a too low flux and should therefore be relatively 
unimportant and free-free absorption instead dominate. In figure we show an 
updated version of the figure in fl^| . The most interesting point is that most 
Type Ib/Ic SNe, SN 1983N, SN 19941 and SN 1998bw, fall into the high velocity 
category, while Type IIL SNe, like SN 1979C and SN 1980K, as well as the Type 
Iln's SN 1978K, SN 1988Z, SN 1998S fall in the free-free group. SN 1987A is 
clearly special with its low mass loss rate, but is most likely dominated by SSA 
Q^,[l4|] . Mrk 297A is an apparent radio supernova in the clumpy, irregular 
galaxy Mrk 297 p2[; its supernova type is not known. 



8 Discussion and conclusions 



Circumstellar interaction of supernovae gives an important window on the nature 
of stars that explode and their evolution leading up to the explosion. Mass 
loss rates for the red supergiant progenitors of Type II supernovae range from 
~ 2 x 1(T 6 M yr _1 for SN 1999em to > 2 x 1(T 4 M yr" 1 for SN 
1979C and SN 1986J fsj]. Evidence for CNO processing has been found in a 
number of supernovae, including SN 1979C g|, SN 1987A [|| and SN 1995N 
[ ftlj . In some cases, the reverse shock appears to be moving into gas that is H 
poor and O rich, e.g., SN 1995N [Q; this relates to the total amount of mass 
loss before the supernova. The complex circumstellar environment of SN 1987A 
has become clear because of it proximity (see Chapter by McCray). For more 
distant supernovae, studies of polarization and spectral line profiles can reveal 
asymmetries, as in SN 1998S {52). 

The evidence on circumstellar interaction is especially useful when it can 
be combined with information from other aspects of the supernovae, such as 
their light curves and stellar environments. For example, from the pre-supernova 
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stellar environment of SN 1999em, Smartt et al. |Q deduced an initial mass of 
12 ± 1 Mq. The supernova was of the plateau type, implying that hydrogen 
envelope was largely intact at the time of the supernova. This is consistent with 
the relatively low rate of mass loss deduced for the supernova progenitor J66[ . 

In addition to information on the evolution of massive stars and their explo- 
sions, circumstellar interaction provides an excellent laboratory for the study of 
shock wave physics. Compared to older supernova remnants, the shock velocities 
are higher and the time evolution gives an additional dimension for study, al- 
though there is little spatial information in most cases. VLBI observations can, 
however, in this respect be extremely valuable, as demonstrated by SN 1993J. An 
object where both the spatial and time dimensions are accessible is SN 1987A, 
which has turned out to be an excellent source for the study of shock waves (see 
Chapter by McCray). 

We are grateful to John Blondin for providing Figure 2 and to Peter Lundqvist 
for comments on the manuscript. This work was supported in part by NASA 
grant NAG5-8232 and by the Swedish Research Council. 
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